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Constraints for the constants of hypothetical Yukawa-type 
corrections to the Newtonian gravitational potential are ob- 
tained from analysis of neutron scattering experiments. Re- 
strictions are obtained for the interaction range between 
10~ 12 and 10~ 7 cm, where Casimir force experiments and 
atomic force microscopy are not sensitive. Experimental lim- 
its are obtained also for non-electromagnetic inverse power 
law neutron-nucleus potentials. Some possibilities are dis- 
cussed to strengthen these constraints. 
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1 Introduction. 
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Hypothetical new long-range interactions have been discussed from different 
points of view. First, developments in quantum field theory have lead to the 
possibility of the existence of number of light and massless elementary parti- 
cles [1]. Exchange of such particles between two bodies can reveal itself as 
an additional interaction potential of Yukawa or power-type and result in a 
deviation of the gravitational force from Newtonian law. On the other hand 
the theoretical models have been developed recently with extra spatial dimen- 
sions [2] leading to an additional forces at short distances. The deviations of 
results of measurements of gravity forces at macroscopic distances from calcu- 
lations based on Newtonian physics can be seen in the experiments of Galileo-, 
Eotvos- or Cavendish-type [3] performed with macro-bodies. At smaller dis- 
tances (10~ 7 — 10~ 2 ) cm the effect of these forces can be observed in measure- 
ments of the Casimir force between closely placed macro-bodies (for review see 
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[4]) or in the atomic force microscopy experiments. At even smaller distances 
such experiments are not sensitive enough, and high precision particle scatter- 
ing experiments may play their role. In view of absence of electric charge the 
experiments with neutrons are more sensitive than with charged particles, elec- 
tromagnetic effects in scattering of neutrons by nuclei are generally known and 
can be accounted for with high precision. As an example of precision scattering 
experiments we present here the analysis of the measured total neutron scat- 
tering cross section by heavy nuclei in the eV-keV energy range, and consider 
some possibility to increase the sensitivity in inferring new interactions from 
the neutron scattering experiments. 



2 Non- Newtonian gravity. 

For the potential of the form 



U(r) = -e- r ' x (1) 
r 

where g = —aGMm, G - the Newtonian gravitational constant, M and m are 
the masses of gravitating bodies (nuclear and neutron), the Born amplitude for 
the scattering by homogeneous gravitating ball (nucleus) with radius R is 

2mg\ 2 
ti\q 2 \ 2 + 1) 

Here the nuclear formfactor is 

F{q) = J^IMQR) ~ qRcoa(qR)], (3) 

where R is nuclear radius. 

Inverse power- (A+l) gravitation potentials 

WO = f O", (4) 

where = —ot^GMm^ tq is arbitrary length and may be taken 1 fm, and 
A=l,2,... appear [2] in models with large extra dimensions for r < R c , (R c is 
the characteristic compactification radius). 

The value of the interaction radius A in the Yukawa-type potential of the 
Eq.(l), the power in the power-law potential of the Eq.(4), the compactifi- 
cation radius R c of extra dimensions and the strengths of these interactions a 



/(?) = t2 :;r, m 
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and ajy expressed here as their ratios to the newtonian interaction, are not pre- 
dicted by theory and are the subject of experimental investigation. Restrictions 
for the Yukawa-type potential down to cm range were obtained previously 
from experiments with macro-bodies. The bounds obtained in these experi- 
ments on the value a increase from ~1 at the millimeter distances to 10 35 at 
the distances ~10~ 8 cm. The summary of these results is shown in the Figure 
1 with corresponding referencies. 

3 Low energy neutron scattering - basic relations. 

The problem of inferring the admixture of the weak long-range potentials from 
neutron cross section data is not trivial. During the years of attempts of the 
measurements the values of neutron-electron and neutron electric polarizability 
scattering amplitudes different approaches were used with different success (for 
review of some of these methods see [20]). The most precise data for these 
constants were obtained from the measurements of the neutron total scatter- 
ing cross section by heavy nuclei, the differential cross section generally are 
measured with much lower precision. The contribution of the neutron-electron 
interaction to the neutron-nucleus scattering cross section is proportional to the 
electric charge of the nucleus, the contribution of the neutron electric polariz- 
ability is proportional to the square of the charge of the nucleus, therefore the 
measurements of scattering are performed for heavy nuclei. The contribution 
of non-Newtonian gravity should be proportional to the mass of the nucleus, 
it means that in this case the heavy nuclei should give the largest effect. On 
the other hand analysis of neutron scattering by heavy nuclei presents the most 
serious difficulties in view of great effect of resonances in neutron scattering. 
For light nuclei ( 4 He, C) the resonance effect and the contribution of electro- 
magnetic interactions are much lower but there is no precision data on neutron 
scattering by these nuclei up to now. 

In most publications the formalism used for the extracting the amplitudes of 
long-range interactions is described too shortly so that many important details 
are omitted. We use here the method close to proposed in [21] and most ap- 
propriate for analysis of cross sections for high enriched samples, in which the 
overwhelming contribution comes from one isotope. The experimental values of 
coherent scattering lengths are not used in this approach. The modification in 
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our approach consists in the following: we do not subtract the resonance contri- 
bution from experimental data before fitting the experimental data but include 
all resonance terms into the expression for the scattering amplitude and fit the 
experimental data to calculated cross section. Such procedure seems more reli- 
able because resonance contribution to the scattering amplitude interferes with 
potential neutron-nucleus scattering amplitude and with the amplitudes of the 
long-range interactions. 
Scattering amplitude 



where fo(q) is the amplitude of nuclear s-scattering, fi(q) is the amplitude of 
nuclear ^-scattering, / Sch is the amplitude of electromagnetic (Schwinger) scat- 
tering, f ne is the amplitude of n — e scattering, f p is the amplitude of neutron- 
nucleus scattering due to neutron electric polarizability, f w is the amplitude of 
hypothetical additional long-range neutron-nucleus interaction. 
Transmission cross-section 



°{k) = J \f(q)\ 2 dn « 47r|/ (fc) + f ne (k) + f p (k) + f w {k)\ 2 + 127r|/i(fc)| 2 + 



In this expression f^k) = / f^dQ/An, where "i" stands for "0" , "ne" , >" , 
and "u>"; \f\{k)\ 2 = / \fi(q)\ 2 dVl/l2ii] instead of the terms J\fi\ 2 dQ are used the 
terms | / fidQ\ 2 . The latter is valid for the amplitude fo which does not depend 
on the scattering angle, and for the amplitudes f n - e , f p , and f w which are from 
three to four orders of magnitude lower than nuclear scattering amplitude fo. 
The error induced by such an approximation as was shown by direct calculation 
does not exceed the tenth of millibarn that is significantly below the precision 
of the total neutron cross section by heavy neuclei, which in the best case is 
not better than several millibarns. a a is the capture cross section, a Sch is the 
Schwinger scattering cross section due to the interaction of neutron magnetic 
moment with nuclear electric charge, Aa sol is the correction to the scattering 
cross section due to solid state structure and dynamics effects, and a inc is the 
incoherent cross section [22]. 

We do not neglect here small terms in the scattering cross section due to 
interference of ^-scattering with the amplitudes f p , ((Jip) and f w (ai w ). 



/(?) = fo(q) + fi(q) + f,M + fne(q) + fp(q) + /»(<?) 



(5) 



+CT lp + G\ w + a a + CT Sch + (J inc + A(J sol . 



(6) 



fo 




(7) 
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k = 2.1968 x lO -4 ^^ 1 / 2 is the neutron wave vector (fm _1 ) in center of mass 
system (neutron energy in the lab. system in eV). 
S-matrix 

where summation is performed over known s-resonances, T n j = T n j(Ej)k/kj, 
Tj = Tnj H - r^yj, and Qj are the neutron widths, the total widths, and the 
statistical weights of j th s-resonance. 
Upon introducing 

j kj 4(E — Ej) 2 + r 2 ' { ) 

and 

v =v— g^Vj do) 

^( 2 ) Y% 4 (^-^) 2 + r f 

where r° • are the reduced neutron widths, the real and the imaginary parts of 
the nuclear scattering amplitude are 

Re/o = " cos(2 ^ o) E « - sin ( 2(5 o) E (2 ) (ii) 

and 

Im/o = + cos(2<5 ) £ (2) -sin(2£ ) . (12) 

For the distant and unknown resonances, when E <C Ej>, and r 7J / <C T nj v <c 
Ej/, we use the expansion 

E (i0 « E (i0 (fc = o) + E' (i0 (fc = 0)*; + E^)(fe = o)fc 2 A (13) 

(i'=l,2) so that: 



nj 

and 



_ x ^ yj'L nj'*- jf , \ - yj' 1 nj' , , \- gj^ nj' 1 7/ .2 /ir\ 
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The contribution of these distant resonance terms to the cross section term 
proportional to k is: 

87r(i?E'(i (o) - R 2 E ( *)(o) + E ( i (o) E ao (o) + E (20 (o) E' (20 (o)), (16) 

and the contribution of these distant resonance terms to the cross section term 
proportional to k s is: 

Y (2i? 4 E (20 (o) - 3i? 2 E' ( ' 20 (o) - 4i? 3 E' (10 (o) + 

+3E'(i (0)E(i (0) +3E'( 20 (0)E( 20 (0)). (17) 

The calculation of these terms with sum over all known [23] resonances for the 
sample with isotopic content 208 Pbo.983 207 Pbo.on 206 Pbo.oo6 [18, 19] yields the 
contribution for /c-proportional term less than 10~ 6 b, and for /c 3 -proportional 
term ~10~ 5 b at the neutron energy about 20 keV. 

As at low energies 6 = —kR' n , where B! n is the nuclear part of the scattering 
radius, expansion of the first term in Re/o gives: 

Re/o = ~ R ' n + 1^ 2 ~ h^ k * + 3i5^ 6 ~ dk^ 8- 

-cos(2<5 ) E(i) -sin(2^ ) E (2) + hk 2 . (18) 

The last two terms represent the contribution from the tails of distant and 
unknown (including at negative energies) resonances, and comes from the ex- 
pansion of the term cos(2^ ) E(i')- At low energies E <C \Ej/\, and Tji <C \Eji\ 
this contribution is: 

gfTnji j2 9j'^nj' g _ ^ [ ^j^2 

ji 2kj'Eji ji 2kj'Ej/ 

The first term enters into unknown scattering radius R', the second term we 
have to consider as unknown (fitted) parameter. 

The amplitude of the neutron-atom scattering due to the neutron-electron 
interaction is 

f ne = -b ne Z(f-h), (20) 
where b ne is neutron-electron scattering length, 



/(*) = -(\/(l + *) - 1), x = 12(fc/</ ) 2 , qoiim 1 ) = 1.9 x W^Z 1 ' 3 

(21) 
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is the atomic formfactor, and 

S(fc) = 1 _(*M + 2 (^_... (22) 

is the nuclear charge fromfactor [24]. 

R N = 1.2027^ 1/3 = 7.127 fm for 208 Pb is the nuclear electric radius. 
The neutron electric polarizability scattering amplitude 

f P = f P o- Q, (23) 

where 

f P o = ®n& 2 ^- = 32.78 mfm (24) 

tl Hn 

for Z=82 and the neutron electric polarizability a n = 1.2 x 10~ 3 fm 3 , the form- 
factor 

Q = \~ l( kR N) + ^(kR N f - —(kR N ) A ... (25) 
The ^-scattering amplitude is 

h = ^(1 - Si)> (26) 

s ^^-^ iE-Si*A (27) 

Here summation is performed over known p-resonances, Tj — ^nj ~^~~^'yji cmd 
Tj, T n j, T 7 j, and gj are the total width, the neutron width, the gamma-width 
and the statistical weight of j th ^-resonance. 



T n j — — T n j(Ej)— , (28) 

where 



(k 7 R) 2 (kR) 2 

^ = TTTW' Uy = ITPF' (29) 



The phase of ^-scattering 

{ k (kR) 2 
■31 + (kR) 



S 1 = -kR + arctg(fci?) + arcsin[- i ; /7 m2 (i? - #!)], (30) 



where is the channel radius, R[ is the p-wave scattering radius. 
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Interference terms of the ^-scattering amplitude and the polarization ampli- 
tude are: 

a lp = 4 7 rRe/ 1 (/c)/ p (-^ - |c 2 ), (31) 

where x = kR. 



Obtaining of constraints for hypothetical long-range 
interactions from neutron scattering experiment. 



The method to search for the admixture of long-range interaction in the neutron- 
nucleus scattering was proposed by Thaler [25] in application to the measure- 
ment of the neutron electric polarizability. The idea consists in the search for 
non-even in the neutron wave vector terms in the power series of the experi- 
mental angular distribution or the total neutron cross sections in keV energy 
range. 

Indeed, in result of interference of the Born scattering amplitudes for the 
long-range potentials of the Eq.(4) and the short range nuclear scattering am- 
plitude, the following terms appear in the low energy wave vector k expansion 
for the total cross section (Eq.(6)): 



flw = 

1 



2k 



for N = 1, 
for N = 3, 
for N = 5, 



f*w = -^a 2 br 2 MkR) foriV = 2, 
h w = 7^a A br 4 k 2 ln(kR) for N = 4, 
f 6w = -^a 6 b4ln{kR)k A for N = 6. (32) 



hw = -^asbiirlk 
fow = -^a^brlk 3 

Here b = b = 1.3 x 1(T 35 fin" 1 for the nucleus 208 Pb. 

The potentials for gravitating balls with radius R for r> R may be obtained 
by integration of the potentials of the Eq. (4) over the ball volume and have the 
form: 



N=l, U 1 (r) = 



2R 3 r 



R 2 - r 2 .R + r. ' 
m|— 1 + Rr 



N = 2, U 2 (r) = 



2 r 



AT = 3, U 3 (r) = 



2i? 3 



R\ 

3r ! 



T 

-ml 
2 l R 



R-r 1 
R + r , 



R 



2r Ul R-r l 



R 



R 2 



(33) 
(34) 
(35) 
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* = W = ,, ( ,,f r " fl2)2 ; (36) 

5 

^ = 5 > ^(r)= (r2 ^ 2)3 ; (37) 
7\r a tt ( \ 96Tq R 2 + 5r 2 

N = 6, C/eM = 5r (r 2_ i?2)4 - ( 38 ) 



Born scattering amplitudes for these potentials 

R+e 



In = / VN(r)sm(qr)rdr, (39) 
9 J 



where Vat = Un/qn, e is the cut-off parameter in extra-gravity interaction, it 
is taken at the level of electroweak scale ~10~ 4 fm [2]. These amplitudes have 
very complicated view, we are interested in the low neutron energy expansions 
of these amplitudes. 

For example, if to omit in these expansions the constant and even in the wave 
vector terms (which is impossible to distinquish from contribution from short- 
range nuclear interaction) we have the following contributions to the formfactor 
in the total scattering cross section: 

U w = \a A brl\\n(kR) + ln(fce)] • k 2 (40) 

D 

for the potential of Eq.(36), 

f 5w = ^a 5 brl[l + ^(kR) 2 ]-k" (41) 
15 7 

for the potential of Eq.(37) and 

few = -^a 6 br b [ln{kR) + In(fce)] • k 4 (42) 

for the potential of Eq.(38). The first terms in the brackets correspond to the 
pointlike potentials of Eqs.(4), the remaining terms come from the fact that 
the sources of potentials are spherical bodies. For the potentials of Eqs. (33-35) 
with N<3 no additional terms appear in these expansions for the scattering 
amplitudes for gravitating balls in comparison to the cases of the pointlike 
scatterers of Eqs. (32). 
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As an example of the measurement of the total neutron cross section in eV- 
keV energy range we used two sets of data from the publications [18, 19] and 
[26]. These groups measured the total cross section of highly enriched 208 Pb 
samples with atomic concentations: 208 Pbo.983 207 Pbo.on 206 Pbo.oo6 in the energy 
range from 1.26 eV up to 24 keV with declared precision from ~ 20 mb to 2.5 
mb in [18, 19] and 208 Pb .97i2 207 Pbo.oi79 206 Pbo.oi08 in [26] with similar precision. 

The outlined formalism was used for the procedure of extracting the con- 
tribution of possible new long-range interactions to the total neutron cross 
section. The MINUIT program [27] was used with three fitted parameters: 
the s-scattering radius R' n , the term corresponding to the contribution to the 
s-scattering amplitude from distant and unknown resonances hi, and contribu- 
tion from the unknown long-range potential f w . The constraints obtained from 
both sets of data: [18, 19] and [26] were very close. 

The neutron-electron scattering length was varied between b ne = —1.32 x 
10~ 3 fm and b ne = —1.59 x 10~ 3 fm [28-32], the value of the neutron electric 
polarizability was varied in the limits ±30% around a n = (1.22 ± 0.19) x 10~ 3 
(fm) 3 , the latter value being obtained as an averaged over the results from two 
recent measurements of the Compton scattering by deuteron [33, 34]. 

Figure 1 shows obtained limits for non-Newtonian gravity from these neutron 
scattering data in terms of relative potential strength a in Eq.(l) as a function 
of the Yukawa length scale A. It shows also the present status of constraints 
from the previous larger scales experiments. 

Figure 2 shows the limits on the value of the inverse power law potential of 
Eq.(4), parametried in the form: 

U(r) = Uo(lfm)M(—) N+ \ (43) 

r 

where M is the relative nuclear mass. 

It is hardly believable that radical increase in precision of keV-neutron scat- 
tering cross section may be achieved in the forthcoming years. Analysis of sys- 
tematic effects and uncertainties in the processing data of neutron scattering 
arising mostly from resonance structure of cross section shows that the limiting 
precision lies in the region 5 mb. Not the least is the problem of correct account- 
ing the background in the time of flight measurements [18, 19, 26, 35]: different 
methods of measurement and subtracting the background yield significantly 
differing results. 
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Possibly some increase in sensitivity in the search for new interactions in the 
neutron-nuclei scattering may be achieved by measuring the energy dependence 
of the asymmetry of slow neutron scattering by noble gases - the method of 
Fermi-Marshall [36, 30]. 

The neutron-atom scattering amplitude 

a(0) = a N + a n - e Zf{6) (44) 

in result of interference of the nuclear and electron scattering leads to the scat- 
tering cross section; 

= a 2 N + 2a N a n - e Zf(9), (45) 



dtt 

with the relative value of the anisotropic term 

all <2/v 

Atomic formfactor may be calculated with sufficient precision according to [24] 

/(9) = [l + 3(«-)T 1/2 , (47) 

with 

q = pZ 1 / 3 (A- 1 ). (48) 

In these expressions is the nuclear scattering amplitude, a n _ e is the ampli- 
tude of n-e scattering, Z is the charge number. The value of (3 is taken from 
the tables [37]. 

Significant diffraction effect appears in the neutron scattering by noble gases 
due to atom-atom interactions [38-43]. The static structure factor depends on 
interatomic interaction potential U (r) and atomic density n: 

oo 

S(q) = 1 + n / (e- u ^ kT - l)e^ r dh. (49) 
o 

For the spherically symmetric potential U (r) the structure factor is: 

A 00 

S(q) = 1 + n— J ( e - u{r)/kT - l)sm(qr)rdr. (50) 
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The scattering asymmetry is introduced as the ratio of intensities of the scat- 
tered neutrons to the angles 9\ and 62, (usually 9\ + 62 = n). 

5 <*> (51) 



Figure 3 shows the results of the calculations of the effects of hypothetical in- 
teractions on the center-of-mass scattering asymmetry by Xe gas together with 
the effects of neutron diffraction and neutron-electron interaction. The curves 
show the results of calculation of asymmetry for the forward and backward 
angles, respectively, 30° and 150°. For Xe several best approximations for the 
potential of interatomic interaction were taken from the literature. 

5 Constraints for hypothetical long-range interactions 
from antiprotonic atoms. 

Some constraints on long-range forces may be obtained from the measurements 
of energies of electric-dipole transitions in hadronic atoms. The average distance 
of hadron from nucleus is 

< r >=i — 2 — (fm) > (52) 

where n and / are the quantum numbers of the state, and m r is the reduced 
mass of the bound hadron. This distance is in the range of dozens to thousand 
of fermies for studied antiprotonic atoms. The most convenient transitions 
between the states (n,l) are of the type (n, n — 1) — > (n — l,n — 2) with n 
large enough to decrease as much as possible the effects of strong interaction 
on the energy states of orbiting hadrons. To estimate the contribution to the 
energy shifts of these states from the potentials of Eq.(l) and Eq.(4) we have 
to calculate the expectation value in the hydrogenlike atomic states of these 
potentials: They are 

AE n = - 2 — ( ) , 53) 

where Z is nuclear charge, cio = h 2 /fie 2 , for the potential of Eq. (1) and 

AEn = ^ ( 2gro(. +1) (2n- j V-l)! 
r nao (2n)! 
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for potentials of Eq.(4). The difference of experimental energies of transitions 
E(n, I) — E(n — 1, / — 1) and the calculation of the QED predictions for these 
transition energies may be assumed to be the effect of long-range interaction. 
No real difference was observed in the most precise experiments we used for 
obtaining these restrictions. We used the reported values of experimental un- 
certainty 4xl0~ 7 eV for transitions (32,31)^(31,30) in antiprotonic helium [17] 
and the value 50 eV for transitions (5,4)— >(4,3) in antiprotonic sulfur [16] as the 
boundaries for the difference in the energy shifts due to these long-range inter- 
actions. Figure 1 shows these constraints as a function of A for the potential of 
Eq.(l). For the potential of Eq.(4) similar constraints may be obtained assum- 
ing that the compactification radius Rcomp ^>< r > - the average distance 
of the orbiting antiproton from the nucleus. The latter is calculated according 
to Eq.(51) and is ~ 0.2 A for the case of antiprotonic helium [17], and ~ 50 
fm for the case of antiprotonic sulfur [16]. The constraints for the potential 
of Eq.(4) are q?2 > 6xl0 33 for the measurement [16] and ct^ > 6xl0 32 for the 
measurement [17]. 

The possibility of probing additional space dimensions with fast neutron 
scattering experiments was considered recently in [44]. I am grateful to Dr. W. 
M. Snow (Indiana University) for attracting my attention to this paper. The 
author is grateful also to Dr. A. B. Popov for important comment and to Dr. 
G. S. Samosvat for consultations and discussions. 
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Fig. 1. Bounds olq vs. A (90% C.L.) from different non-Newtonian grav- 
ity experiments. The area above the lines is excluded: 1 - Cavendish- 
type experiment of Seattle group [5]; 2 - microcantilever Cavendish- 
type experiment of Stanford [6]; 3 - Casimir force measurements; con- 
straints from[7, 4], based on earlier experiments of [8, 9]; 4 ~ Casimir force 
measurement [10], the bound from [11]; 5- Casimir force measurements 
with atomic force microscope[12], bounds from [4]; 6- neutron quantum 
levels in the Earth gravitational field, experiment [13], bounds from [11]; 7 
- van der Waals force measurements [9], bounds from[14, 4]; 8 - atomic- 
force microscopy experiment [15]; 9 - bounds obtained in this work from 
X-ray energy measurements in antiprotonic sulfur atoms[16]; 10- bounds 
obtained in this work from the measurements of the transition frequencies 
between large n levels in antiprotonic He-atoms[17]; 11 - bounds obtained 
in this work from analisys of neutron total cross section scattering by 208 Pb 
nuclei[18, 19]; 

Fig. 2. Limits (90% C.L.) on the value of the interaction potential Uo(l/ra) 
vs N in the expression U(r)=U (lfm)M(^r + : 1 - from X -rays m 
p— 4 He - atoms; 2 - from X-rays in p—S - atoms; 3 - from the measure- 
ment of the total n— 208 Pb - scattering; 4 - the same as 3, but additional 
interference terms for scatterting amplitudes for gravitating balls of Eqs. 
(40-42) were included to the expression for the total scattering cross sec- 
tion of Eq.(6). 

Fig. 3. The calculated neutron energy dependence of the center-of-mass 
scattering asymmetry (6\ = 71/6,62 = 57r/6) from 1 atm pressure Xe gas 

target in result of different effects: 1 and 2 - n e interaction with the 

values of n e scattering lengths -1.32-10~ 3 fm and -1.59-10 _3 fm; 3, 4, 5, 

and 6 - diffraction effect for different approximations of the interatomic 
Xe-Xe interaction: 3 and 4 - two approximations of Lennard- Jones poten- 
tial for Xe[38, 39], 5 and 6 - potentials of more complicated form from [40]; 
7, 8, 9 and 10 - effect of hypothetical potential of the form Eq.(l) for A 
equal to 10 5 fm, 10 4 fm, 10 3 fm and 10 2 fm and a equal to 10 23 , 10 25 , 10 27 
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and 10 9 respectively; 11, 12 and 13 - effect of hypothetical potential of 
the form Eq.(4) for N = 1 and compactification radius 10 5 fm, 10 4 fm and 
10 4 fm and a\ equal 10 29 , 10 30 and 10 31 respectively. 
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